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One of the outstanding problems in J/U physics is a systematic understanding of 
the differential photo-production cross section dcrldz{'y -|- p —> J/tp X), where 
2 : = in the proton rest frame. The theoretical prediction based on the non- 

relativistic QCD (NRQCD) factorization formalism has a color-octet contribution 
which grows rapidly in the endpoint region, 2 : —> 1, spoiling perturbation theory. In 
addition there are subleading operators which are enhanced by powers of 1/(1 — z) 
and they must be resummed to all orders. Here an update of a systematic analysis 
is presented. The approach used to organize the endpoint behavior of the photo¬ 
production cross section is based on a combination of NRQCD and soft collinear 
effective theory. While a final result is not yet available, an intermediate result 
indicates that better agreement between theory and data will be achieved in this 
framework. 
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1. Introduction 

The production of J/i/' in high energy electron-proton interactions has 
proven to be a rich subject with great potential for furthering our un¬ 
derstanding of the perturbative and non-perturbative regimes of QCD. In 
particular a number of different kinematic regions are accessible in ep col¬ 
lisions so that a range of different J/ijj production phenomena can be ex¬ 
plored with a single experiment. The HI and Zeus collaborations have 
gathered a large amount of data on J/ip production and study a diverse set 
of topics including processes where the scattered electron is in the forward 
region so that the photon exchanged with the proton is real. 

The theoretical framework within which the production of non- 
relativistic bound states can be systematically treated is the NRQCD fac¬ 
torization formalism This approach relies on a non-relativistic effective 
theory of QCD which is a systematic expansion of QCD about the 
limit of vanishing quark velocity: u ^ 0 . The resulting effective field the¬ 
ory consists of an infinite number of operators each scaling as a definite 
power of u, however, at at a given order in v there are only a finite number 
of terms. As a consequence working to a specified numerical accuracy only 
requires including those contributions that affect the calculation at the level 
of significance desired. 

Reference [ 1 ] postulates that the J/'tj) photo-production cross section in 
ep collisions is of the form 

du{'^-\-p J/V’ + AT) = fi/pi^dd{i+j cc[n]+X)i^F{cc[n] J/ip) , ( 1 ) 

where dd{i -I- 7 ^ cc[n] -I- A) is the short-distance scattering cross sec¬ 
tion for producing a cc pair with non-relativistic relative momentum in a 
color and angular momentum state indexed by n. The parton distribution 
function fi/p gives the probability for finding parton i in the proton, and 
F{cc[n] —>■ Jftp) gives the probability for the cc pair with quantum numbers 
n to hadronize into a J/ip- There is an implied sum over quantum numbers 
n, with each F{cc[n] J/ip) scaling as with 7(71) the scaling dimen¬ 

sion of the operator. The short-distance coefficient dd{i -I- 7 ^ cc[n] -I- X) 
scales with powers of the strong coupling constant as, and can be enhanced 
by kinematic factors as well. As a result there is a competition between 
suppression by and enhancements of dd, and extra care must be taken 
when deciding which terms to keep in Eq. ( 1 ). In particular there are two 
effects that can greatly enhance the short-distance cross section: fragmen¬ 
tation contributions and color-octet contributions ®. 

In Refs. [ 7 , 8, 9 ] the total J/ip photo-production cross section (7(7 -|- 
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p ^ J/^jj + X) is calculated and compared to HERA data as well as other 
experimental data. The conclusion is that theory and data agree as long 
as the non-perturbative NRQCD production matrix elements are allowed 
to be negative. However, given the interpretation of the NRQCD matrix 
elements as the probability for a cc pair to fragment into & J/ip and any 
number of soft particles there is a prejudice against negative values of the 
matrix elements. Furthermore there is uncertainty whether the NRQCD 
factorization formalism is valid for the total cross section since the J/'ll) is 
not necessarily produced with large transverse momentum. References [7, 
9] also considered the differential cross section da/dz{^+p J/'tj){z) +X) 
where 0 = Pp • Ptf; /Pp • p-y for J/ip transverse momentum greater than 1 GeV. 
The results of Ref. [7] compared to HI data are shown in Fig. 1. The 



Z 


Figure 1. The J/ 1/1 energy distribution dujdz at the photon-proton center-of-mass en¬ 
ergy yfs = 100 GeV integrated in the range py > 1 GeV from Ref. [7]. 


theoretical curve is an increasing function of z while the experimental points 
indicate a flat distribution, and it is clear that the color-octet contribution 
dominates the differential cross section at larger values of z. The rise in the 
cross section is due to terms of the form 


(^) 


1 -z 


+ 


( 2 ) 


which arise at next-to-leading order in the color-octet differential cross sec¬ 
tion. These plus-distributions become arbitrarily large as z ^ 1, and cause 
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a breakdown of perturbation theory. As a consequence such terms must be 
resummed to make a sensible comparison of theory to data in the large z 
region. 

In addition to the breakdown of perturbation theory, the non- 
perturbative NRQCD expansion breaks down when z ^ 1, because of a 
set of subleading operators which scale as — z)". This was first ex¬ 

plored in Refs. [11, 12], where the infinite tower of enhanced operators was 
resummed into a shape function. In Ref. [13] an analysis of the differential 
cross section da/dz including the shape function was carried out. The spec¬ 
trum compared to HI and Zeus data is shown in Fig. 2. The introduction 



z 


Figure 2. The J/ip energy spectrum at y/s = 100 GeV integrated in the range pj_ > 1 
GeV compared to HERA data. Solid, dashed-dotted, and dashed lines are the color- 
singlet contribution plus the color-octet shape function contribution for three different 
values of a shape function parameter. The dotted line is the color-singlet contribution 
alone. Plot from Ref. [13]. 


of a shape function tames the endpoint divergences of the color-octet con¬ 
tribution at NLO, however, agreement with data is not great. In particular 
the shoulder in the theoretical prediction is too large and peaked too far to 
the right. This is due to the inconsistent treatment of the plus-distribution 
terms which also need to be resummed before a comparison of theory and 
data can be made. In Ref. [15] it was shown that in e+e” ^ J/ip + X at the 
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endpoint perturbative resummation in addition to the shape function sig¬ 
nificantly broadens the differential distribution relative to either including 
only the shape function or only carrying out a perturbative resummation 
without shape function. 


2. The breakdown of the NRQCD expansion 


It is helpful to study the kinematics ofy-l-p^ J/il^ + X as z ^ Ito gain a 
better understanding of the physics that leads to the breakdown of NRQCD 
factorization in this kinematic regime. In the proton-photon center-of-mass 
frame 

Pj = Pp = Pcc = + k >^, ( 3 ) 

where = (1,0,0, —I), = (1, 0, 0,1), M = 2mc, and and are the 

4-velocity of the J/ip and the residual momentum of the cc pair in the J/ip 
respectively. In terms of the scaling variable z the J/if) velocity is 

P'i = +P± + ( 4 ) 

where Mj_ = ~ p\^ with By momentum conservation 


Px =^7 


Pp - Pcc 

^ lyiM 
2 \M 


xVs 


1 - 


MMl\ 

szM^ J ^ 




Setting fc = 0 the invariant mass squared of the final state is 


Px = 


sM / 


1 - 


MMl 


2 


( 6 ) 


szM^ 

This is much larger than A ~ 1, GeV as long as M^/M — z ~ 1. Note 
the requirement for p\>Q translates into p^ < (sz — MM^){M^/M — z) 
so p± is forced to be small at the endpoint. Thus away from the endpoint 
the decay products can be integrated out and NRQCD factorization holds. 
However if z ^ /M the invariant mass of the decay products becomes 
much smaller than the h • px momentum component and the final state 
consists of a J/ip and a jet in the direction of the incoming hadron. In this 
kinematic regime NRQCD factorization no longer holds, and an effective 
field theory which describes the light-like final state must be used. The 
appropriate theoretical framework is a combination of NRQCD to describe 
the J/ip and soft collinear effective theory (SCET) I6,i7.i8,i9 describe 
the light-like decay products. A combination of SCET and NRQCD has 
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already been used to study a variety of Quarkonium decay and production 
processes in kinematically restricted regions 20 , 21 , 2245 , 23 , 24 , 25 , 26 ^ 

Before diving into a detailed description of SCET it is useful to deter¬ 
mine the region in z where SCET is applicable. To do this return to Eq. (6) 
and include the residual momentum of the cc pair introduced in Eq. (3). 
The scaling of the residual momentum can be found by boosting from 


the J/tp rest frame where ~ 
frame. In light-cone coordinates 

fc'' ~ A 


A to the proton-photon center of mass 


M 


M.i, 


(7) 


where z ~ 1 has been used. Note the plus light-cone component is enhanced 
by while the minus component is suppressed by the same amount. 

Including the residual momentum in Eq. (6) gives 

MMl\ 


2 ^ 

M^\M 


1 - 


, m 2 2 MMl 


n-k + ... (8) 


where suppressed pieces are dropped. The last term in the equation above 
scales as sA/M^ so the final state becomes sensitive to nonperturbative 
physics when z > (M^ — K)/M. 


3. Soft Collinear Effective Theory 


SCET describes the dynamics of highly energetic particles moving close to 
the light-cone interacting with a background field of soft quanta. The in¬ 
teraction of the collinear particle with the background introduces a small 
residual momentum component into the light-like collinear momentum so 
that collinear particles have momentum + k^, where k^ ~ A. 

However, collinear particle do not only interact with the soft background, 
they also couple to other collinear particles. As a consequence the SCET 
Lagrangian consists of two sectors: soft and collinear. The SCET La- 
grangian was first derived in what is called the label formalism 16 , 17 , 18 , 19 ^ 
and was subsequently formulated in position space illustrate some 

important properties of SCET consider the quark sector of the Lagrangian. 
It can be split into two pieces: one coupling collinear to soft 




.in 




(9) 


where P,n,p is the collinear quark field labelled by collinear momentum = 
n-p/2+p'P, n = (1, 0, 0,1), and in- D = in-d + gn-Ag, with As the soft 
gauge field. This expression looks very much like the HQET Lagrangian 
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with the velocity replaced with the light-like vector and was first 
derived in Ref. [29]. The second piece of the collinear Lagrangian consists 
of interactions of only collinear particles among themselves 

•Tc = ■ ^n,q + ) (10) 

where g is the collinear gauge field, and the collinear covariant derivative 
is iD^ = +igAn^q. The operator projects out collinear label momen¬ 
tum: V^^n,p = n-p/2+p^^)^n,p The SCET Lagrangian is invariant 
under separate collinear and soft gauge transformations which provide a 
powerful restriction on the operators allowed in the theory . Furthermore 
the Lagrangian is invariant under a global U{1) helicity spin symmetry, 
and must be invariant under certain types of reparameterizations of the 
collinear sector of the Lagrangian 30,3i 

A remarkable consequence of the gauge symmetries of SCET is the 
factorization of soft and collinear effects. Towards this end the soft Wilson 
line is introduced 

E(a;) = Pexp^fg J ds n ■ Aus{ns)^ , (11) 

and the collinear fields are redefined as follows: 

^n,p = TC® , (12) 

After the field redefinitions the soft gluons decouple from the collinear fields, 
i.e. Ccs ^ 0 in Eq. (9), and the collinear Lagrangian becomes independent 
of soft physics. At higher orders in the SCET expansion this decoupling 
does not occur, and factorization is broken. 

4. Factorization at the Endpoint 

The soft-collinear factorization properties of SCET can be used to obtain 
a factored form for the J/tp production cross section at the kinematic end¬ 
point. An important condition for the derivation to be valid is that a 
sufficient range of p± is integrated over. To be precise we must smear over 
a range p± > AM. If this is not done, or if only a range p_L ~ A is inte¬ 
grated over the intrinsic transverse momentum of the partons in the proton 
must be taken into account 32,33,34 jg g.^ important point because an 

experimental cut of p± > 1, GeV is usually implemented in HERA data. 
This cut lies in the regime where non-perturbative momentum is important 
and may introduce a sensitivity to the transverse momentum of the partons 
in the proton. 
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There are two steps in factorizing the photo-production differential cross 
section at the endpoint. The first step is matching QCD onto SCET where 
the collinear particles have a typical off-shellness of order AM. This for¬ 
mulation of SCET is called SCETi. In the next step SCETi is matched 
onto SCETn where collinear fields have a typical off-shellness 0{h?') To 
be concrete one of the two leading color-octet contributions to J/'0 photo¬ 
production at the endpoint is considered: a cc pair produced in a color-octet 
configuration which hadronizes to J/ip via a chromomagnetic spin-flip 
transition. 

The derivation is easiest in the parton-photon CM frame where 

K = = xpp = + Pp + > (^3) 

with the J/ip approximately at rest, = (1,0,0, 0). First the QCD current 

is matched onto the SCETi current 


LO 

J^ix) = (x), (14) 


where u) = ujnP and P = Vn^f2 + V±. At lowest order in as{M) 






(15) 


The second step is to integrate out the scale AM by match the differ¬ 
ential cross section in SCETi onto SCETn. To derive the cross section it is 
best to take a step back and consider the differential cross section in QCD: 


drr — C 

= J dSe-^P--y(p\JlmJ/i^ + X){J/i^ + X\My)\p) 

(16) 

Inserting the expression in Eq. (14) gives the SCETi cross section 


^ ^1,2 

X {p\[xlpSf^^Bli;pm\J/i; + X){J/^P + X\\iPl5f^^Blx-p\{y)\p) , (17) 


and soft-collinear factorization properties of SCET outlined in the previous 
section can be used to obtain a factored form for this cross section: 


da 

dz 


aasC^ 

ttMzs 


{OtCSo)) / dn-kS 


1 - 


M 

M,i, 


z +n-k J(jj^{n-k). (18) 
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Here the momentum space functions and S are Fourier transforms of 
collinear and soft matrix elements respectively: 

/ r/r? • k 

(19) 

and 

S^ico) = 2M{OtCSo))S^{to) 

= I ^e-i-^-y{0\xlpYT^Y^p{0)V^i;l,YT^Y^X-p'{n-ym.{20) 
By construction f dujS(oj) = 1. 

Next SCETi is matched onto SCETn- This entails performing an op¬ 
erator product expansion (OPE) on Ju}^{n-k), where the large scale is set 
by the invariant mass of the collinear degrees of freedom in SCETp 

J^{n-k) ^ 

x{pu\ Tr{r^H”(0)}4p+...+ Tr{T^Hi^(0)} \pu), (21) 

where corrections are suppressed by 0{A/M). The coefficient function is 
dimensionless, and therefore must be a function of the ratio of the invariant 
mass squared and the large momentum component squared. Furthermore, 
at leading order in as{AM) the perpendicular momentum of the J/'ip is 
0{A) so that the labels opj: must be zero. The SCETn operator above 
can be related to the familiar parton distribution function which gives the 
probability of finding a gluon in the proton 

i^(ml Tr{T«i?f(0)}4p+...yTr{T^i?f(0)} \pn) 

spin ^ 

= -n-pp dxS{u;+-2xn-Pp)fg/p{x), (22) 
Jo 

where fg/p (x) is the parton distribution function. Using this result the final 
factored form of the cross section is 

J=M|C'^“(M)|"^4^)^ I d^S40Cn{^-l + z-A)fgJ^y 

( 23 ) 



February 2, 2008 4:51 Proceedings Trim Size: 9in x 6in 


SFleming 


10 

5. Summing Logarithms 

Large logarithms are summed using the renormalization group equations 
(RGBs), which is a two step procedure in this case. First the effective 
theory currents in Eq. (14) are run from the scale fxn = '/s to the scale 
= M\/l — 2 using the SCETj RGEs. Next the SGETn operators in 
Eq. (23) are run to the scale /x ~ A. Details of the running are left for a 
later publication Here preliminary results are presented in Fig. 3, which 
shows the differential cross section (in units of the total color-octet ^S'o 
photo-production cross section). The solid line includes both perturbative 
resummation and non-perturbative resummation i.e. the shape function. 
The dashed line includes only perturbative resummation. This should be 
compared to Fig. 2 where only the shape function is included. In the case 
where either the shape function alone is used, or only perturbation theory is 
resummed the spectrum is too sharply peaked to be compatible with data. 
However, when both the shape function and perturbative resummation is 
included the spectrum softens considerably, which gives hope that theory 
will be compatible with data in a complete analysis. 
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